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Abstract —In order to foster renewable energy integration, 
improve power quality and reliability, and reduce hydrocarbon 
emissions, there is a strong need to deploy energy storage systems 
(ESSs), which can provide a control medium for peak hour 
utility operations. ESSs are especially desirable at the residential 
level, as this sector has the most untapped demand response 
potential. However, considering their high acquisition, operation, 
and maintenance costs, isolated deployment of ESSs is not eco¬ 
nomically viable. Hence, this paper proposes a sharing-based ESS 
architecture, in which the demand of each customer is modeled 
stochastically and the aggregate demand is accommodated by 
a combination of power drawn from the grid and the storage 
unit when the demand exceeds grid capacity. The optimal size 
of energy storage systems is analyzed an analytical method is 
developed for a group customers with single type of appliances. 
This framework is also extended to any network size with 
arbitrary number of customers and appliance types, where the 
analytical method provides a tractable solution to the ESS sizing 
problem. Finally, a detailed cost-benefit analysis is provided, 
where the results indicate that sharing-based ESSs are practical 
and yield significant savings in terms of ESS size. 


I. Introduction 

Over the past few years, the power grids have become 
more stressed due to the steady increase in peak demand and 
more fragile due to the integration of intermittent renewable 
energy resources. Such a negative trajectory poses a profound 
threat to power grids, leads to a rise in hydrocarbon emissions, 
and exposes new economic challenges. One effective solution 
to mitigate the aforementioned issues, is to deploy energy 
storage systems (ESSs), which can act as an energy buffer 
and decouple the time of generation and demand by storing 
cheaper and cleaner off-peak hour electricity and delivering it 
during the peak load periods m. Integration of such storage 
units has multi-faceted monetary benefits for different players. 
For instance, end-users can enjoy the reduced electricity cost 
and lessened interruptions, while utility operators can enjoy 
improved grid reliability and asset utilization, and possibly 
deferral of system upgrades. Moreover, policymakers benefit 
from reduced hydrocarbon emissions and the mainstream 
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integration of renewables, which are closely linked to energy 
security m and 0- 

While employing storage devices at microgrids has certain 
benefits, their deployment based on the current technologies 
may still not be economically viable 0. Therefore, optimal 
sizing of the storage units based on the realistic needs of the 
grids is a critical step towards the efficient operation of the 
grid. Specifically, over-provisioning ESS size entails costly 
and underutilized assets, whereas under-provisioning reduces 
its operating lifetime (e.g., frequently exceeding the allowable 
depth of charge level degrades its health). Hence, there is a 
strong need to develop analytical models to solve the sizing 
problem. Furthermore, the sizing requirements for storage 
units are usually application-dependent. For instance, energy 
arbitrage applications for independent system operators require 
100s mega watts (MW) in size, while storage units employed 
for energy management for communities or microgrids require 
10s kilo watts (kW). In this study, our main focus is on the 
latter one, as the residential sector has a great potential for 
peak demand reduction in the U.S. 0. 

ESS sizing has received some attention in the literature. The 
work in a presents a sizing approach for single industrial 
customers for peak saving applications. The sizing problem 
is solved by maximizing the net benefit, which is the sum of 
reductions in the electricity bills minus the operation costs and 
the one-time acquisition cost. Similarly, a proposes a sizing 
framework using similar cost models for a microgrid, but it 
also considers savings due to the storage of energy generated 
from renewable resources. The work in 0 develops a sizing 
approach based on stochastic network calculus to size the 
storage units coupled with solar power generation. It employs 
loss of load probability as the main performance metric to 
provision the resources. Another probabilistic approach is 
presented in 0, which couples forecasting errors in wind 
generation with the storage unit and sizes the storage unit 
according to a desired level of accuracy. Moreover, 0 solves 
the sizing problem by using stochastic dynamic programming 
to minimize the operation cost and employs the storage device 
for balancing the load in a wind farm application. From 
the power engineering point of view, the sizing problem is 
usually solved via simulations m-m. However, simulations 
techniques are usually computationally expensive, and their 
accuracy depends upon the availability of data traces. 

In this paper, we develop an analytical framework for 
optimal energy storage sizing for sharing-based end-user ap¬ 
plications. The proposed framework contributes to the existing 
literature in the following ways: 

« The existing analytical methods for storage sizing focus 
on settings with one customer, or are developed for 
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renewable integration and do not explore the underlying 
user dynamics. The proposed framework can cope with 
any network with arbitrary number of consumers with 
different levels of demand. Our analytical results show 
that a sharing-based design of storage units exhibits 
substantial gains over user-level designs. 

• The existing methods for multi-consumer settings are 
simulation-based. Lack of closed form expressions for 
the storage size-cost model limits the applicability of the 
model in optimizing the operations. The advantage of the 
proposed analytical method is hat it establishes the exact 
optimal sizing, and subsequently, is computationally less 
expensive. 

• Majority of the existing literature assumes that residen¬ 
tial users employ storage units for energy management 
applications without considering the associated costs. 
However, the studies conducted by Electric Power Re¬ 
search Institute Cl and Sandia Laboratories El reveal 
that the cost of employing stand-alone ESS is very higlfl 
Hence, we provide a detailed economic analysis for ESS 
applications in residential usage and show that storage 
technologies can be economically viable, only if operated 
in a shared manner. 

Proper sizing of storage units is ultimately linked to the cus¬ 
tomer electricity consumption model. Several measurement- 
based studies m-m, and m show that electricity con¬ 
sumption of households can be represented by On-Off models. 
For instance, the work presented in m measures the elec¬ 
tricity consumption of 20 households at fine-grained intervals 
(every 6 second) and models the loads with continuous¬ 
time Markov chains. Driven by this observation, we adopt a 
Markovian fluid model to represent consumer’s demands, and 
the analyzes rely on the stochastic theory of fluid dynamics. 
We establish the interplay among the minimum amount of 
storage size, the grid capacity, number of consumers, and the 
stochastic guarantees on outage events. 

We note that studying storage units at a network level is 
of paramount significance as they are expected to become 
integral to smart energy grids. More specifically, ESSs will 
be employed at smart residential, business complexes, and 
university campuses, to name a few, to reduce peak hour 
consumption. Clearly, in such sharing-based applications, the 
size of the energy storage is linked to the customer population, 
the power drawn from the grid, and the load profile. This 
relation is highly non-linear due to multiplexing gains which 
are computed by the percentage of reduction in the required 
amount of resources with respect to the baseline case of 
assigning peak demand to each user. 

II. System Description 

We consider a community of consumers in which the 
demands of N users are accommodated by the power grid 
capacity along with a shared energy storage system unit, 
in which capacity fluctuates over time. The grid capacity at 

The cost/benefit regime depends on different assumptions such as avail- 
ability of solar rooftops, energy consumption statistics, project duration, and 
utility tariffs, to name a few. See references for details 


TABLE I: Notations 


Parameter 

Description 

Ct 

Power drawn from grid at time t. 

N 

Number of users in the single class case. 

N 

Vector of the number of users for K classes. Note that this is not the 
number of houses since in one house there can be multiple appliances 
requesting demand. 

Rk 

Demand of a user of class k. 

Afc 

Arrival rate of charge request for class k. 

Mfc 

Mean service rate for the customer demand of class k. 

B 

Size of the energy storage unit. 

S(t) 

ESS depletion level, 0<.S(t)<B. 

L n (t) 

Aggregated load on the system when n users are On. For single class 
index i is used. 

Fi (x) 

Steady state cumulative probability distribution function of ESS 
charge level. 


Grid power allocated per source (C/N). 

K 

ESS per user ( BIN). 

Wfc 

Effective demand for customer class k. 

C 

Bustiness parameter. 

V 

Round-trip efficiency of the storage unit. 


time t is denoted by Ct and ESS depletion level is S(t), for 
t £ R + . Furthermore, the energy storage unit has the following 
parameters: 

« Energy rating or the size of the storage is denoted by B 
(kWh). 

• Power rating is the rate at which storage can be charged 
or discharged. The charging rate is assumed to be P c < 
C t . Vt, and discharge rating Pd is related to the energy 
rating B since B P,i x (desired support duration). 
Desired support duration is typically equivalent to the 
length of the peak hour. 

• The efficiency of a charge-discharge cycle is modeled by 
r] £ [0,1] to capture the percentage of stored energy or 
the fraction of energy that is transferred after losses are 
extracted. 

• Dissipation losses represent a small percentage of losses 
that occur due to leakage. For simplicity in notations, 
dissipation losses are ignored. 

• For economic analysis, we follow the standard assump¬ 
tions made in JT] and assume that the project will be sited 
for 15 years, with 10% discount rate and the economic 
analysis is carried out by computing the net present value 
of the storage in the first year by using the cost per kW 
or kWh. 

Furthermore, since the consumers do not necessarily have 
identical demands, we consider K customer classes, which 
are distinguished by the amount of electricity demands, rep¬ 
resented by where Ilk denotes the energy demand per 

time unit for customer type k. Let Nk denote the number of 
consumers of type k and the vector N represents the number 
of consumers of each type, that is N = (Ni, ■ • •, Nk). As 
established in urn ini, m. uni and fT31 the consumption 

pattern of each consumer can be well-represented by a two- 
state “On/Off” process. We define the binary variable s l t k to 
represent the state of consumer i of type k at time t such that 

ik ( 1 consumer i is On 

— [ 0 consumer i is Off ' * 

When a customer is in the “On” state, it initiates an energy de¬ 
mand. The duration of demand is modeled statistically, which 
is adopted to capture the variety types of consumers’ demands. 
Specifically, the duration of type k customer’s demand is 
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assumed to be exponentially distributed with parameter pk- 
Furthermore, we assume that the requests, that are transitions 
from “Off” to “On”, are generated randomly and according 
to a Poisson process with parameter A&. Hence, for each 
consumer i of type k at any time t we have 


p (4* = i) 


Afc + pk 


( 2 ) 


In order to formalize the dynamics of ESS, we define Lik(t) 
as the charge request of consumer i of type k. Then, there are 
exactly three cases that define the rate of change in the storage 
unit: (i) the storage can be in the fully charged state, and the 
aggregate demand is less than grid power C t \ (ii) the storage 
can be completely discharged, and the aggregate demand is 
more than the grid power C t ; and (iii) any partially charged 
state with any customer demand. Therefore, for the rate of 
change in the storage level of the ESS, we have 


dS(t) 

dt 


0, 


0, 


v(Ct 


t)), 


if S(t) = B & 


K N k 


E E vLik(t) < c t 

k=l i= 1 


if S{t) = 0 & 


K N k 


E E vLik{t) > c t 

k =1 i= 1 


otherwise 


(3) 


Due to stochasticities involved (in consumption and gener¬ 
ation), by choosing any storage capacity B, only stochastic 
guarantees can be provided for the reliability of the system 
and always there exists a chance of outage, which occurs 
when available resources fall below the aggregate demands 
by the consumers. By noting that 5(f) denotes the energy 
level that the storage unit needs to feed into the grid to avoid 
outage, an outage event occurs when the necessary load from 
the storage unit exceeds the maximum available B. Hence, 
we define e-outage storage capacity, denoted by B(e), as the 
smallest choice of B corresponding to which the probability 
of outage does not exceed e G (0,1), i.e., 


% J min B ... 

B “ { s.t. P (St >B)<£- (4) 

Our goal is to determine the e-outage storage capacity B{e) 
based on grid capacity Ct , the number of users TV, and 
their associated consumption dynamics. For the simplicity of 
mathematical expressions, we scale the storage parameters and 
instead of B/r) we redefine B as the maximum amount of 
energy that can be stored, and similarly Pd and P c represent 
the actual power ratings. The notations are summarized in 
Table U We start our analysis for deriving B{e) for a single 
class case. 



Fig. 1: Composite model for N independent users for single user 
type (k=l). Each user becomes active (“On”) at rate A and becomes 
inactive (“Off”) at rate p. The aggregate demand depends on the 
active number of users. 

Since we have N independent consumers each with a two- 
state model, by taking into account their underlying arrival 
and consumption processes, the composite model counting 
the number of users accessing the storage unit at a given 
time can be modeled as a continuous -time birth-death process. 
Specifically, this process consists of (N + 1) states, in which 
state n G {0,..., N } models n consumers being active and 
accessing the storage unit, i.e., 

N 

state at time f is n if E = n , (5) 

i=1 

and drawing nR p units of power from the storage unit. As 
depicted in Fig [I] the transition rate from state n to state n +1 
is (N — n )A and, conversely, the transition from state n + 1 
to state n is (n+ l)p. Hence, for the associated infinitesimal 
generator matrix M, in which the row elements sum to zero, 
for i, n G {0,..., N} we have 


M[i, n] 


-((N - i )A + ip) 
ip 

(N - i) A 
0 


n = i 

n = i — 1 & i>0 
n = i + 1 & i < N 

otherwise 


(6) 


By denoting the stationary probabilities of state n G 
{0,..., N} by 7 r„ and according defining 7r = [7To, tt\, ..., 7Tjv], 
these stationary probability values satisfy irM = 0. 


B. Analyzing Distributions 

Given the dynamics of accessing the storage unit, in the next 
step we analyze the statistical behavior of the ESS charge level. 
Specifically, we define x) as the cumulative distribution 
function (cdf) of the ESS charge level when i G {0,..., N} 
consumers are depleting the storage unit at time t, i.e., 

N 

Fi(t,x ) = P(5(t) < x and s{ = i) . (7) 

i =i 

Accordingly, we define the vector of cdfs as 


III. Storage Capacity Analysis For Single Class 
Customers (K = 1) 

A. Storage Access Dynamics 

When the grid can serve all the consumers’ demands, there 
will be no consumer served by the storage unit. On the 
other hand, when the grid capacity falls below the aggregate 
demand, the consumers access the storage unit. Since the 
requests of the consumers arrive randomly, the number of 
consumers accessing the unit also varies randomly. 


F{t,x) = [F 0 (t,x) , Fi(f,x) , ... , F N (t,x)] . (8) 


Based on this definition, the next lemma delineates a differ¬ 
ential equation which admits the cdf vector as its solution and 
is instrumental for analyzing the probability of outage events, 

i.e., P(£ ” = L i(t) > C t + B). 

Lemma 1. The cdf vector Fit. x) satisfies 


dF(t, x) 


■ D = F(t, x) ■ M 


dx 


(9) 
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where D is a diagonal matrix defined as 

D A diag [-Ctn , (1 - C t )n , ... , (N - C t )p\ , (10) 

and matrix M is defined in ©. 

Proof: 

In order to compute the probability density functions, we 
find the expansion of Fi(t,x) for an incremental change At 
in t, i.e, Fi(t + At,x). Note that during incremental time At, 
three elementary events can occur: 

1) one inactive consumer might become active, i.e., i 
increases to i + 1; 

2) one active consumer becomes inactive, i.e., i reduces to 
* — 1; or 

3) the number of active consumers remains unchanged. 
Since the durations of arrival and departure of consumers are 
exponentially distributed as cdf Fi(t,x) can be expanded 

Fi(t + A t,x) 


= [TV - (i - 1)] ■ (AA t) • Fi_i (t , x) 


one consumer added 

+ [i + 1] • (ftAf) • F i+1 (t,x) 

'-v-' 

one consumer removed 

+ [1 — ((TV — i) A + ip) Af] • Fi(t, x — (i — Ct) • pAt) 

S * v ^ 

no change 

+ o (At 2 ) , (11) 

where o (At 2 ) represents the probabilities of the compound 
events and tends to zero more rapidly than Af 2 (and Af) as 
Af —> 0. Next, by passing the limit 

Fi(t + At, a;) 


lim 

At —^0 


At 


it can be readily verified that (fill simplifies to 

dFi(x,t) 


dt 


= [N — (i — 1)] - (A) - Fi_i(t,x) 
+ [i + 1] • (ft) • F i+1 (t,x) 

- [(TV - i)X + ip] ■ Fi(t, x) 

^ ^ dFi(t,x ) 

- (« - C«). (M)-^— , 


( 12 ) 


where we have defined F_i(f,x) = Fjy +l (f, x) = 0. By 

recalling that the design is intended gor a long-term steady- 
state operation we have dFi{x,t) / dt = 0. Hence, (IT2l > can be 
rewritten as 

r)F (t X ^ 

(i - Ct) • (m) ■ = [N-(i- 1)] • (A) • F^x) 


dx 


+ [« + !]■ (ft) • Fi + \{t, x) 


- [(N-i)\ + in]-Fi(t,x) . (13) 

By concatenating all the equations ( fl3] > for all i £ {0,..., N} 
we obtain the compact form 


dF(t, x) 
dx 


• D = F(t, x)-M . 


(14) 


The solution of the first order differential equation given 
in (fl4l) can be expressed as a sum of exponential terms. The 
general solution requires computing (TV + 1) eigenvalues of 
the matrix MD~ X and the general solution can be expressed 
as US): 

N 

F(t, x) = ^2,ai <t>i e-xjp{zix) , (15) 

2=0 

where z.j is the i th eigenvalue of MD~ X with the associated 
eigenvector </>,, which satisfy Zif>iD — <piM. The coefficients 
{ao, ..., ctjv} are determined by the boundary conditions, e.g., 
Fi(t, 0) = 0 and F,(f, oo) = 1. 

In order to compute the probability distribution in (fl~5l >. we 
need to determine the eigenvalues of MD~ X , the eigenvectors 
and coefficients {at}. Note that, since x > 0 and 
Fj (f, x) is upper bounded by 1, all of the positive eigenvalues 
and the corresponding a, must be set to zero, hence this 
reduces the computational complexity and (031 ) simplifies to: 

F(t,x)= ^2 a i <t>i exp(zjx) . (16) 

i:Re[zi< 0] 

It can be further observed that since Zi<f>^D = (f)jM, one of 
the eigenvalues must be zero. Then by setting zo = 0, the cor¬ 
responding eigenvector can be computed from 4>qM = 0. We 
also earlier showed that the steady state probability distribution 
7r of the TV + 1 state Markov chain can be computed from the 
same equation, that is 7 xM = 0. Since, the eigenvector cf o is 
known and one of the eigenvalues is zo = 0, we can write 
4 >o = 7r. Therefore, (IT6l > further simplifies to fl9l : 

F(t,x) = n+ ^2 a i <t>i exp(zix) . (17) 

i:Re[zi< 0] 


C. Single User Storage Capacity (K = 1, TV = 1) 

For computing the desired e-outage storage capacity B(e) 
by leveraging the cdf vector found in ([T71 > we start by a simple 
network with one user (TV = 1). The insights gained can be 
leveraged to generalize the approach for networks with any 
arbitrary size TV. When TV = 1 the infinitesimal generator 
matrix M defined in © is 


M = 


(18) 


-A A 

. M -A 1 J 

For finding the expansion of F(t, x) as given in ( fl~5T > we need 
to find the eigenvalues of MD~ X , i.e. zo and zi, where D is 
defined in (ITOt . Based on (ITOt we find that 


MD~ l = 


i 

c t 


l-C, 


1 -Ct 


Hence, the eigenvalues are 
zq = 0 and 


X 


1 -Ct 


(19) 


( 20 ) 


where we have defined \ — “• It can be readily verified that 
the eigenvector associated with zi is 0i = [1 — Ct , Ct}. 
Therefore, according to (fTTI) we have 


F(t,x) = tv + ai 0i exp(zix) . 


( 21 ) 
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Finally, by finding the coefficient a± we can fully characterize 
F(t,x). This can be facilitated by leveraging the boundary 
condition F\(t, 0) = 0, which yields 


Fi(t, 0) = m + cnCt = 0 , (22) 


where we have that iri 


-r^—. Therefore 


Oil 


X 

Ct(l+x) ’ 


(23) 


which, subsequently, fully characterizes both cdfs Fo(t, x) and 
Fi(t,x) according to 


F 0 (t,x ) = 7r 0 + «i(l - C i )exp( 2 ;ia;), 

and Fi(t,x) = tt\ + a±Ct exp(zi;c) . 

As a result, by recalling the definition of Fi(t,x) in CD*, the 
probability that the storage level St falls below a target level 
x is given by 


P(5 t < x) = F 0 (x) + Fi(x) = 1 + a\ exp(zix) . (24) 


Given this closed-form characterization for the distribution of 
St, we can now evaluate the probability term 

P (S t > B) , (25) 


which is the core constraint in the storage sizing problem for¬ 
malized in (Q}. Specifically, for any instantaneous realization 
of Ct denoted by c we have 


contributors to the probability distribution, ||20l shows that, 
the asymptotic expression for F\r(x) is given by 

1 / u 

Fn{x)= 2 V */(?)(?+ A(1-?))JV (27) 

x exp (-Nifd) - gd) x ) 
x exp(-2 v / {fd)d + A(1 - s))Nx}), 


where, 

/(c) - log 

a ?(1 + A) — A 
W= cd-A) * 


c(l + A) — A 
? + A(l-?) ’ 


¥>(<0 = ? log(?) + (1 - ?) log(l - ?) - f log(0 + log(l + A) 

■0(1 -?) 


ad) = z + 0.5 (? + A(1 - 0) ■ 
A(l-2?) 


m 


z a (1 _a) + . 


and ip = 


(c + A(l-c)) ’ 

(2^-l)(c(l + A)-A) 3 
? (l-c) 2 (? + A(l-s)) 3 ' 


In this set of equations, time is measured in units of a single 
average “On” time (1 ///). Furthermore, k and <r are defined 
as the ESS per user ( B/N ) and the grid power allocated per 
one source, respectively. Furthermore, we denote the variable 
v as the power above the mean demand allocated per user as 

W = ^TTX- 


P {S t >B) = 



P (S t >B\C t =c ) f Ct (c ) dc 

oli exp(ziB) fc t (c) dc 


c t 



X 

c( 1 + X) 


exp 




f Ct {c) dc . 


Therefore, by noting that Zi = % — j s negative, the 
probability term P(S't > B) becomes strictly decreasing in 
B. Hence, the smallest storage capacity B that satisfies the 
stochastic guarantee l! :, (,S) > B) < e has a unique solution 
corresponding to which this constraint holds with equality. In 
the simplest settings in which grid capacity Ct is constant c 
we find 


B{e) 


c(l - c) _ log ec( 1 + x) 

X-Xc-c x 


(26) 


D. Multiuser Storage Capacity (K = 1, N > 1) 

In this subsection we provide a closed-form for the prob¬ 
ability term P(5't < x) for arbitrary values of N , which we 
denote by F^(x). Computing all F^(x) terms through com¬ 
puting their constituent terms Fi(t,x), especially as N grows, 
becomes computationally expensive, and possibly prohibitive 
as it involves computing the eigenvalues and eigenvectors 
of MD~ X . By capitalizing on the observation that for large 
number of users N the largest eigenvalues are the main 


IV. Storage Capacity Analysis For Multi-Class 
Customers (K > 1, N > 1) 

In this section, we consider the case where the storage serves 
more than one customer type. We start our analysis by noting 
that the continuous-time Markov chain model presented in 
Fig. □ becomes a K dimensional Markov process. Recall 
that total number of users of each class is represented by 
N = (Ni, ■ ■ ■ , Nk) and we let n(t) = \n\it)... 
represent the number of sources of type k that are On at 
time t. Furthermore, n(t) represents the state space of the 
Markov process and similar to single class case given in ©, 
the transitions can occur between neighboring states, and the 
transition matrix M = {M(f,n)} is given by, 

M ( n , A +(n)) = ( N k - n k ) A fc , (28) 

M(n,A~(n))=n k p k , (29) 

M ( n, n') = 0, otherwise. (30) 

where 

A+(m, • • • ,n k , ■ ■ ■ ,n K ) = (ni, ■ ■ ■ ,n k + 1, • • • ,n K ) 

A”(m, • • •, n fc , • • •, n K ) = (ni, • • •, n k - 1, ■ • •, n K )• 

In the single class case, the transitions can only occur be¬ 
tween 2 neighboring states. In the multi-class case, however, 
the transitions can occur in 2 x k different neighbors. For 
instance, assume that there are two classes, k = 2, and the 
Markov chain is in state (1,1) (i.e., one active user from 
each class). Therefore, there exists four possible transitions 
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(1,0), (0,1), (2,1) and (2, 2), and the transition rates depend 
on the (Afc,/Zfc) pairs. To that end, the transition rates given in 
( l28l > represent customer arrivals, whereas the transitions rates 
in ( [29b account for new customer arrivals, and similar to single 
class case there can not be any transitions beyond neighboring 
states. This is given in (l30l) . 


A. Effective Demand 

The analysis in Section IIII-BI provided for the single¬ 
class settings is also valid for the multi-class case. However, 
computing the eigenvalues becomes infeasible when K > 1. 
Hence, we follow the decomposition method proposed in ITT) 
to efficiently evaluate the eigenvalues z t and the coefficients 
a.i. Our goal in this section is to compute u> k which is a 
deterministic quantity acting as a surrogate for the actual 
aggregate stochastic demand l22l . In other words, assigning 
u>k amount of resources to each customer class will satisfy ©. 

We define parameter £ = 1 — e x ! B £ [0,1] to account for 
the susceptibility of the storage to customer demand pattern. 
This parameter will be an essential part of u> k . First we 
will consider large storage units, then we will show that the 
proposed scheme is also valid for any storage size. Assume 
that we are able to choose N = (N±, • • •, Nk) in a way that 
ensures P(5* > B) < e for small e and the eigenvalues are 
relabeled such that zq > z\ > ... > 0, then the following 
theorem holds. 


Theorem 2. Let = {N : P (St > B) < e}. For large 

storage size B and for small e, we have, 


lim 

B=oo,e —>-0 


logg 

B 


—■► C € [ oo, 0] 


Furthermore, let 


and, 


B=\N:J2 0J k(ON k <c\, 


B=\n :J2M0Nk<C l 


where effective demand can be computed by EB, 


fRk + fife + A/t - y (C Rk + Bk — Afc) 4- AXkPk 

^(0 =- 57 - 

(31) 

Then, B C B{B, e) C B. 

Proof: For a given set of users, N, the following holds. 


p (St > B) __ J2 k a k(l T (f k )exp(z k B) 

£ £ 


where z k ’s are the non-positive eigenvalues of MD~ 1 , <fi k 
are the corresponding eigenvectors, and l 7 = [1,1, - - -, 1]. 
Recall that, the coefficients a k can be obtained by solving the 
following set of equations. 


J2k a kfik = 0, 

= 0 , 

aol^o = 1, 


if D k > 0, 
if z k > 0, 
otherwise. 


TABLE II: Cost of poor power quality. HD 


$/kW 

Average 

High 

Residential 

$0.10 

$0.60 

Small C&I 

$0.42 

$2.52 

Large C&I 

$1.42 

$14.00 


TABLE III: Average outage cost [25l 


$/kW 

15 Min. 

30 Min. 

1 Hour 

2 Hours 

Residential 

$0.05 

$0.60 

$2.60 

$3.95 

Small C&I 

$8.65 

$16.01 

$23.37 

$48.91 

Large C&I 

$4.79 

$7.46 

$10.12 

$17.96 


In fZB, it is shown that Zq > z k , for \/k > 1, and Zq satisfies 
C = 'Bk a k{zo)N k . Then the following holds, 

— —- = a o l T 0oexp (2:o_c)B (l + o(l)), (32) 

£ 

as B —>• oo. Differentiation shows that, effective demandu>k(y) 
decreases as y increases. Thus, for N £ B, Zq < f and 
ipf s >b') 

v * ~ —- -A 0 as B ^ oo, which indicates that N £ B(B,£). 
In a similar manner, when N (£ B, then zo > f and 

jp| O ) 

— - -A oo, and therefore N ^ B(B,e). This completes 
the proof. ■ 

The interpretation of Theorem[2]is that © holds if and only 
if (jj k (f )N k < C. Furthermore, the following remarks hold. 

k 

Remark 2.1. The corresponding effective demand for Q = 
0 is u>k{0 = x^+ii k ’ is the mean customer demand 

for each class. Another interpretation of this result is that, 
before an outage event occurs at the storage, each customer 
turns “On” and “Off” their appliances many times that the 
aggregated demand equals to the mean customer demand. This 
result serves as a lower bound for the effective demand. 

Remark 2.2. At the other boundary, the corresponding effec¬ 
tive demand for f = z is its peak demand R k . In this case, 
each user is at “On ” at all times, hence peak demand must 
be allocated. 

Remarks given above suggest that the system is more 
susceptible to bustiness (longer “On” durations) for larger f 
values. Also, similar to the multi-class case, the largest con¬ 
tributor to the probability distribution is the largest eigenvalue, 
i.e., Zq. 


V. Economies of ESS Deployment 

The proposed framework targets enhancing energy manage¬ 
ment via deploying for end users and the monetary associated 
benefits of ESS deployment ©. 

1) Improved power quality refers to voltage sags and out¬ 
ages experienced by customers, which in many cases 
remains unnoticed. However if occur for sufficiently 
large durations, they damage customers’ appliances. The 
costs of momentary outages, according to a survey study 
conducted in 1 241 . are summarized in Table HI1 

2) Improved power reliability refers to the usage of storage 
units during outages and blackouts. In order to quantify 
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TABLE IV: Average Time Of Use (TOU) rates in the United States 

m 



Summer 

Winter 

$/kWh 

Peak 

Off-Peak 

Peak 

Off-Peak 

Residential 

$0.25 

$0.06 

$0.13 

$0.06 

Small C&I 

$0.18 

$0.05 

$0.12 

$0.05 

Large C&I 

$0.06 

$0.04 

$0.05 

$0.04 

TABLE V: Demand Charges 


Summer 

Winter 

$/kW-Month 

Peak 

Off-Peak 

Peak 

Off-Peak 

Residential 

$0.00 

$0.00 

$0.00 

$0.00 

Small C&I 

$15.00 

$15.00 

$8.00 

$8.00 

Large C&I 

$12.00 

$12.00 

$10.00 

$10.00 


the cost of power reliability, we adopt the measure 
used in (25), i.e., the cost per kW shown in Table [III] 
Obviously the impact of the outages depends on the 
duration and the frequency of the events. We adopt the 
statistics conducted by Q] and assume that the consumer 
average interruption duration in one year is 88 minutes. 

3) Reduced Time of Use (TOU) charges include the savings 
occurred by eliminating the use of peak hour electricity, 
and using storage units. “Time Of Use” tariffs may vary 
in different territories and for different seasons. In Table 
IIV1 we present the average TOU tariffs in the U.S. 

4) Demand charges is a significant portion of the commer¬ 
cial and industrial customers’ bill. Some utilities also 
apply these charges to residential customers. It is usually 
computed by the amount and the duration of the peak 
usage a. Since it is not easy to compute these charges 
without knowing the load profile, we did not consider 
this in our calculations. An overview of demand charges 
is presented in Table [V] 

It is noteworthy that our goal in discussing this example sce¬ 
nario is to show how our framework can make the ESS deploy¬ 
ment economically viable. Using ESSs is mostly application- 
dependent and there can more benefits if storage is coupled 
with distributed generation options. Additional benefits will 
include profits by energy trading and reduced electricity 
bills. Also, the cost of carbon emission decreases (due to 
the elimination of peaking generators). The exact amount of 
savings depends on the specific policies and regulations. For 
commercial and industrial customers, there will be additional 
benefits due to demand charges. 

VI. Numerical Examples 
A. Single Class Customers 

In this section, we provide numerical examples to explain 
the system dynamics and show how the proposed framework 
can be used in typical peak shaving applications. In the first 
group of setting, uniform customer demand is considered 
(K = 1). We use the aforementioned normalized values (unit 
time is measured in //A 1 and unit demand is measured in peak 
demand - R p ). We start by exploring the relationships between 
the number of users ( N ), ESS size (in f? p /r _1 units) and 


the corresponding underflow probability for a given system 
capacity, which remains constant over time, C. Charge request 
rate per single user A is set to 0.3, and the mean capacity above 
the mean demand per user is set to v = 0.035 (assumptions 
are valid for all single class case studies). Therefore, the total 
system capacity becomes C = 0.2658IV units. In Fig. |2(a)| 
sizing problem is evaluated for user population N ranging 
from 100 to 350. Considering the fact that, one household 
can employ 10 — 12 appliances, this interval is chosen to 
represent a multi-dwelling building or a typical building on a 
university campus. These findings can be used in the following 
ways. First, system operators can provision the storage units 
and provide statistical guarantees (underflow probability) to 
their customers for a given user population N. For example, 
for a large-scale electric vehicle charging lot (e.g., located in 
shopping mall or airports l26l l with N = 350 charging slots 
in order to accommodate 99.95% of the customer demand 
the ESS size should be selected as B = 7 x R p x /r _1 = 
7 x 6.6 x 2 = 92.4 kWh (considering level-II chargers and 30 
minutes as the unit time). It is noteworthy that the required 
storage size decreases as the user population increases. 

Next, we consider the case where the system operator 
employs an already acquired ESS of sizes B = lOkWh, 
15kWh, or 20kWh. In this case, the critical step is to calculate 
the minimum amount of power to draw from the grid so 
that performance guarantees can be achieved. This case is 
evaluated in in Fig. |2(b)| for N = 200 appliances. For instance, 
suppose ESS size B = lOkWh is already acquired and 
the goal is to meet 95% of the demand, and then system 
operator should draw 0.26 x 200 = 52kW from the grid. The 
primary motivation for the employment of the ESS is to reduce 
the stress on the grid and improve the utilization of power 
system components in a cost effective manner. Hence, our 
last evaluation is on computing the amount of savings in ESS 
size with respect to current common sizing practice, allocating 
peak demand. The results are depicted in Fig. |2(c)| show that 
instead of sizing the ESS to meet the entire customer demand, 
just by rejecting a few percentage of customers, considerable 
savings in the storage size can be achieved. 

B. Multi-Class Customers 

Next, we consider the case where customers can request 
different demand levels, hence different customer classes. In 
this subsection, our main objective is to compute the Uk 
(effective demand ) parameter for each class such that Uk 
replaces Rk and the overflow probability targets are met. We 
start by discussing a toy example. Suppose that users can 
generate four different levels of demands, with the following 
parameters A = {0.3,0.5,0.7,0.9} , fi = {1,1,1,1}, and 
R k = 0.2,0.4,0.6,0.8 (in kW). Then for the target outage e = 
1CT 4 , and ESS size B = lOkWh, the corresponding effective 
demand becomes uj = {0.0515, 0.1567, 0.2958, 0.4550}. 
Alternatively, the system operator can find the minimum stor¬ 
age size with respect to the available grid resources and target 
performance metric. For the simplicity of presentation, the 
remaining numerical results consider two classes of customers 
with the following parameters Ai = 0.5, /ri = 1, If = 0.5kW 

















(a) Community storage sizing for different 
user population. 



(b) Community storage sizing for varying 
grid power, N = 200 appliances (users). 



(c) Savings in ESS Size for varying user 
population 



50 100 150 ,200 

Nji # of Acceptable Class-1 Cust.(R =0.6kW) 


(a) Admission Region for K = 2. 


Fig. 2: Evaluation of single-class case. 



(b) Effective Demand Calculation for K = 2. 
Fig. 3: Evaluation of multi-class case. 



(c) ESS Sizing for K = 2. 


and A2 = 0.7, // 2 = 1, i ? 2 = lkW. Our first evaluation shows 
the admission set for C = 50kW and B = lOkWh, for a target 
overflow probability of £ = 0.001. Results presented in Fig. 
|3(a)| shows that station operator can choose from any set of 
customer numbers from the shaded region. Obviously, since 
the demand for Class 1 customers is less than the other class, 
the system operator can accept more customers from Class 1. 
As a second evaluation, we calculate the Effective Demand for 
varying storage size. This time the target underflow probability 
is set as £ = 0.0005 and the customer population is chosen 
as Ni = 100 and N 2 = 45. The results depicted in Fig. |3(b)| 
shows that instead of provisioning the system according to 
peak demand (Ri = 0 .6kW, li-> = lkW), the use of effective 
demands (e.g., for B = lOkWh, u>\ = 0.25 and w 2 = 0.5) 
reduces the provisioning of resources tremendously. Our final 
results are on the ESS sizing for multi classes. For the same 
set of parameters, we evaluate the storage size with respect to 
power drawn from the grid. Similar to the single class case, 
this result can be used to size the storage unit for a given grid 
power, or it can be used to compute the required grid resources 
for a given storage size. 

C. ESS Economic Analysis 

In this subsection, we provide several numerical examples 
with actual real-world scenarios to show how the proposed 
framework reduces the cost of the customers. In our eval¬ 
uations, we choose the average cost per residential customer. 
We assume that peak appliance demand is set to R\ = 1.2kW, 


each user employs 15 appliances, and there are 250 days of 
peak usage in one year. For simplicity, the peak hour duration 
is assumed to 1 hour. All the cost/benefit calculations are 
normalized to peak hour demand. We consider three cases and 
compare the cost of system operation for different customer 
population and the frequency of appliance usage. 

• Case 1 - Grid Power Only: In this setting, all of the 
peak hour demand is met by the power demand. This 
case reflects the current state of affairs. The cost includes 
the cost of peak charging ($0.25kWh) and all of the 
listed benefits except the demand charges. Obviously, 
since there is no storage unit, users suffer from power 
quality and reliability, and high peak-hour charges. 

• Case 2 - ESS Only: In this case, we consider the case 
where all of the peak hour demand is met by the storage 
unit. The cost includes the cost of off-peak hour charging 
and the storage cost. In this case, all of the listed benefits 
except the demand charges contribute to reducing the total 
cost of the system. 

• Case 3 - Proposed Framework: This is the case where 
we compute the storage size based on the mathematical 
proposed framework. We set the probability target as 
£ = 0.001. By noting that, there can be more than one 
{C,B) combinations that satisfy the overflow targets for 
the evaluations, we assume that grid power is 20% more 
than the average customer demand. 

According to Electric Power Research Institute [|jj, the greatest 
cost of such of energy management and Time-Of-Use pricing 
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(b) ESS becomes viable when more than 11 
customers share the cost for A = 0.45. 
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(c) ESS becomes viable when more than 7 
customers share the cost for A = 0.65. 


Fig. 4: Economical evaluation and cost analysis of the proposed framework in US Dollars ($) per user per one month. 


applications stems from the acquisition and the Operation & 
Management costs of the storage unitfl Generally, the lifetime 
of a storage unit is assumed to be 15 years with 10% discount 
ratf] In this paper, we assume that the unit cost of the storage 
unit per year is $1500, which is the typical cost for a Li-ion 
battery. 

We compute the aggregate cost of each case per user per 
month for three different demand frequencies A = 0.25, 
A = 0.45, and A = 0.65. The results depicted in Fig.[4]have the 
following interpretations. Since the cost of ESS technologies is 
high, the stand-alone ESS setting (case 2) is extremely costly 
and not useful practically, which is in line with the conclusion 
of EPRI m and SANDIA [|2) studies. On the other hand, 
users can enjoy the aforementioned benefits by sharing the 
storage (e.g., in a neighborhood or on a campus etc.), as 
more users share the same resource pool, the cost per user 
decreases. Obviously, the population of users who share the 
same ESS depends on the physical deployments. For instance, 
100 apartments in highly dense regions can be fed by the 
same storage, while this may not be possible in a wide area 
due to losses. Further, as the storage usage increases (higher A) 
the economic benefits improve faster. These evaluations show 
that, it may be practical to share the same storage unit so that 
the total cost of ownership reduces, and customers enjoy a 
variety of benefits. One final note is that the same analysis can 
be applied to the multi-class case by computing the effective 
demands and the corresponding storage size. However, our 
main goal is to conduct the economic analysis per user based. 

VII. Conclusion 

In this paper, we have developed a stochastic analytical 
framework to provision a sharing-based energy storage units, 
which are expected to be employed in peak hour energy 
management systems for residential customers. The analysis 
establishes the interplay among dynamic grid capacity, the 
number of consumers, different appliances types, and the 
guarantee levels for avoiding outage events. We have provided 
a detailed economical analysis and have shown that ESS at 

2 Another highly relevant report is presented in (2j. however we consider 
jil] as it is more recent. 

3 This is used to compute the Net Present Value of the ESS. See [I] for 
further details. 


residential level is economically beneficial if employed in a 
sharing-based architecture. 
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